The Bologna Open Cluster Chemical Evolution (BOCCE) Project: 
midterm results from the photometric sample 
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ABSTRACT 

We describe a long term project aimed at deriving information on the chemical evolution of 
the Galactic disk from a large sample of open clusters. The main property of this project is that 
all clusters are analyzed in a homogeneous way, to guarantee the robustness of the ranking in 
age, distance and metallicity. Special emphasis is devoted to the evolution of the earliest phases 
of the Galactic disk evolution, where clusters have superior reliability with respect to other 
types of evolution indicators. The project is twofold: on the one hand we derive age, distance 
and reddening (and indicative metallicity) interpreting deep and accurate photometric data with 
stellar evolution models, and, on the other hand, we derive the chemical abundances from high- 
resolution spectroscopy. Here we describe our overall goals and approaches, and report on the 
mid-term project status of the photometric part, with 16 clusters already studied, covering an 
age interval from 0.1 to 6 Gyr and Galactocentric distances from 6.6 to 21 kpc. The importance 
of quantifying the theoretical uncertainties by deriving the cluster parameters with various sets 
of stellar models is emphasized. Stellar evolution models assuming overshooting from convective 
regions appear to better reproduce the photometric properties of the cluster stars. 

The examined clusters show a clear metallicity dependence on the Galactocentric distance 
and no dependence on age. The tight relation between cluster age and magnitude difference 
between the main sequence turn off and the red clump is confirmed. 



Subject headings: Galaxy: evolution 
Hertzsprung-Russell (HR) diagram 

1. Introduction 



Galaxy: disk - Open clusters and associations: general 



In the last decade, our understanding of the 
chemical evolution of the Galaxy has significantly 
improved, thanks to the efforts and the achieve- 
ments both on the observational and on the theo- 
retical sides. Numerical chemical evolution models 
nowadays can satisfactorily reproduce the major 
observed features of the Milky Way (e.g. Lacey 
& FaU 1985, Tosi 1988, Matteucci & Frangois 
1989, Ferrini et al. 1994, Giovagnoh & Tosi 
1995, Timmes, Woosley & Weaver 1995, Chiap- 
pini, Matteucci & Gratton 1997, Boissier & Prant- 
zos 1999). They are considered successful only 
when they reproduce the current distribution with 
Galactocentric distance of the SFR (e.g. as com- 



piled by Lacey & Fall 1985), the current distribu- 
tion with Galactocentric distance of the gas and 
star densities (see e.g. Tosi 1996, Boissier & Prant- 
zos 1999 and references therein); the current dis- 
tribution with Galactocentric distance of element 
abundances as derived from HII regions and from 
B-stars (e.g. Shaver et al. 1983, Smartt & RoUe- 
ston 1997); the distribution with Galactocentric 
distance of element abundances at slightly older 
epochs, as derived from PNe II (e.g. Pasquali 
& Perinotto 1993, Maciel & Chiappini 1994, Ma- 
ciel & Koppen 1994, Maciel, Costa & Uchida 
al. 2003); the age- metallicity relation (AMR) not 
only in the solar neighborhood (Twarog 1980) but 
also at other distances from the Galactic center 
(e.g. Edvardsson et al. 1993); the metallicity 
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distribution of G-dwarfs in the solar neighbor- 
hood (e.g. Rocha-Pinto & Maciel 1996); the local 
Present-Day-Mass-Function (PDMF, e.g. Scalo 
1986, Kroupa, Tout & Gilmorc 1993); the relative 
abundance ratios (e.g. [0/Fe] vs [Fe/H]) in disk 
and halo stars (e.g. Barbuy 1988, Edvardsson et 
al. 1993). 

There are, however, several open questions 
which still need to be answered. One of the impor- 
tant longstanding questions concerns the evolution 
of the chemical abundance gradients in the Galac- 
tic disk. The distribution of heavy elements with 
Galactocentric distance, as derived from young 
objects like HII regions or B stars, shows a trend 
in the disk of the Galaxy and of other well stud- 
ied spirals which is usually interpreted as a linear 
negative gradient. It is not clear yet whether the 
gradient slope changes with space and time, and, 
in case, if it flattens or steepens. 

Galactic chemical evolution models do not pro- 
vide a consistent answer: even those that are able 
to reproduce equally well the largest set of ob- 
servational constraints predict different scenarios 
for early epochs. By comparing with each other 
the most representative models of the time, Tosi 
(1996) showed that the predictions on the gradi- 
ent evolution range from a slope initially positive 
which then becomes negative and steepens with 
time, to a slope initially negative which remains 
roughly constant, to a slope initially negative and 
steep which gradually flattens with time (see also 
Tosi 2000 and Chiappini, Mattcucci & Romano 
2001 for updated discussions and references). 

From the observational point of view, the situ- 
ation is not much clearer. Data on field stars are 
inconclusive, due to the large uncertainties affect- 
ing the older, metal poorer ones. Planetary neb- 
ulae (PNc) are better indicators, thanks to their 
brightness. PNe of type II, whose progenitors are 
on average 2-3 Gyr old, provide information on 
the Galactic gradient at that epoch and show gra- 
dients similar to those derived from HII regions. 
However, the precise slope of the radial abundance 
distribution, and therefore its possible variation 
with time is still subject of debate. The PNe data 
of Pasquali & Perinotto (1993), Maciel & Chiap- 
pini (1994) and Maciel & Koppcn (1994) showed 
gradient slopes slightly flatter than those derived 
from HII regions, but Maciel et al. (2003) and Ma- 
ciel, Lago & Costa (2005) have recently inferred. 



from a larger and updated PN data set, a flatten- 
ing of the oxygen gradient slope during the last 
5-9 Gyr. 

We believe that open clusters (OCs) represent 
the best tool to understand whether and how the 
metallicity spatial distribution changes with time, 
since they have formed at all epochs and their 
ages, distances and metallicities are more safely 
derivable than for fleld stars. OCs have gener- 
ally been suggested to show a metallicity gradient 
(e.g., Janes 1979, Panagia & Tosi 1981, Friel 1995, 
Carraro, Ng & Portinari 1998, Friel et al. 2002), 
and Maciel et al. (2005) have recently suggested 
that OC literature data are consistent with those 
on PNe in the indication of a time flattening of 
the gradient. However, the spatial and temporal 
evolutions of the metallicity distribution of OCs 
are still uncertain. 

Bragaglia et al. (2000) using the compilation 

of ages, distances and metallicities listed by Friel 
(1995), and dividing her clusters in four age bins, 
found no significant variation in time of the gra- 
dient slope. This, however, more than reflecting 
the Galaxy evolution, may be an effect of inho- 
mogeneity in the data treatment of clusters taken 
from different literature sources. Carraro et al. 
(1998), using the metallicities of 37 OCs derived 
by Friel's group from low resolution spectra, and 
age and distances based on their own synthetic 
color-magnitude diagrams (CMDs), tried to mea- 
sure the time dependence of the radial metallicity 
gradient. They found hints that the gradient may 
be slightly shallower at the present time than in 
the past, but, again, the results may be affected by 
uncertainties related to the uneven qiiality of the 
literature data they applied the synthetic CMD 
method to. Similarly, Friel et al. (2002) presented 
new, improved data on the metallicity of 39 OCs, 
and derived a smooth radial gradient, with indica- 
tion that its slope may be steeper for old clusters, 
but do not consider this a firm result, due to the 
different space distribution of young and old OCs. 

Against the view of an OCs radial abundance 
gradient, Twarog, Ashman, & Anthony-Twarog 
(1997), and more recently Corder & Twarog 
(2001), argue that what is observed is a step dis- 
tribution in OC metallicity, rather than a smooth 
gradient: clusters with Rqc ^ 10 kpc have solar 
metallicity (with some dispersion), while clusters 
farther away have [Fe/H] ~ —0.3 dex (again with 
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some dispersion). Yong, Carney, & Teixera de 
Almeida (2005) also find a discontinuity at 10 kpc 
in the radial distribution of the OC metallicity, 
and a basement [Fc/H]~-0.5 beyond that radius, 
but agree on the existence of a smooth gradient in 
the inner regions. 

Open clusters are also good tracers of the 
Galaxy formation mechanism (Freeman & Bland- 
Hawthorn 2002), in particular of the thin disk. 
The old OCs may date the Galactic disk, as an 
alternative to White Dwarfs, which still give a 
wide range of values, depending on the adopted 
cooling models (Hansen & Liebert 2003). OGs 
do not offer only a value integrated on the whole 
disk, but give information on its different parts. 
If the discontinuity at 10-12 kpc described above 
is real, it may represent the limit of the original 
thick disk or the connection of OCs to the pres- 
ence of accreted satellites (Frinchaboy et al. 2004, 
2005; Bellazzini et al. 2004) or to the occurrence 
of major merger events (Yong et al. 2005). These 
possibilities are not mutually excluding, since the 
accretion-related clusters would be those in the 
outer part of the Galaxy, where the concept of 
"true" disk is fuzzy. As already noted by Priel 
(1995), old OCs may be the few survivors of the 
original ones that formed in the early Galaxy, or 
may be the product of later accretion or interac- 
tion with merging satellites. 

Finally, OCs are excellent tests for stellar evolu- 
tion theory, since they contain simple stellar pop- 
ulations (i.e. coeval stars with equal metallicity). 
Viable stellar evolution models must therefore re- 
produce the various features (morphology, num- 
ber counts, color and luminosity functions) of their 
CMDs. Comparing theoretical and observational 
CMDs one can discriminate among different mod- 
els and select the most realistic ones. The method 
of using synthetic CMDs to compare with observa- 
tional data has proven to be much more informa- 
tive and rewarding (Aparicio et al. 1990, Tosi et 
al. 1991, Skillman & Gallart 2002) than the clas- 
sical isochrone fitting. These MontcCarlo simula- 
tions allow modeling of several additional parame- 
ters which dictate the distribution of points in the 
CMD, such as stochastic star formation (SF) pro- 
cesses, binary fraction, photometric spread, main- 
sequence (MS) thickness, data incompleteness and 
small number statistics. Consequently, the results 
not only provide a measure of the properties of 



the cluster, but can also constrain the star forma- 
tion history (SFH) and the initial mass function 
(IMF). 

To try to improve our understanding in all these 
subjects it is crucial to build up a large sample 
of open clusters for which the age, distance and 
metallicity are all derived in the most precise and 
homogeneous way, to avoid spurious effects due to 
inhomogeneous treatments. For this reason, we 
arc undertaking a long term project of accurate 
photometry and high-resolution spectroscopy to 
derive homogeneously ages, distances, reddening 
and element abundances in open clusters of var- 
ious ages and Galactic locations, and eventually 
infer from them the evolution of the metallicity 
radial distribution. To this aim wc need to col- 
lect at least 30 clusters, evenly distributed in age, 
metallicity and Galactocentric distance bins to al- 
low for reliable derivation of possible trends (or 
lack thereof). Age, distance, and reddening are 
obtained deriving the CMDs from deep, accurate 
photometry and applying to them the synthetic 
CMD method described by Tosi et al. (1991). 
Accurate chemical abundances are obtained from 
high-resolution spectroscopy, applying to all clus- 
ters the same method of analysis and the same 
metallicity scale. 

Here we describe the results obtained for the 
first half of the photometric sample. Spectroscopic 
abundances have been presented so far only for 
NGC 6819 (Bragagha et al. 2001), NGC 2506, 
NGC 6134, IC 4651 (Carretta et al. 2004), and 
Cr 261 (Carretta et al. 2005) and will be the sub- 
ject of fore-coming papers. This paper is orga- 
nized as follows: the BOCCE photometric sample 
in described in Section 2 and the application of 
the synthetic CMD method in Section 3. The re- 
sults obtainc^d so far by BOCCE on age, distance 
and reddening of the examined clusters are sum- 
marized in Section 4, and discussed in Section 5 
in terms of test for the stellar evolution models. 
Contaminating effects are examined in Section 6. 
A discussion of the mid-term evolution picture ob- 
tained from the currently assembled sample is pro- 
vided in Section 7 and a summary in Section 8. 

2. The photometric sample 

Homogeneity is a fundamental requirement 
when using different data sets to make detailed 
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comparisons between observational data and mod- 
els. To truly understand the variations in the 
different model parameters, one must minimize 
systematic errors. What is needed for the pur- 
poses described above is a large data set of open 
clusters, all observed with the same technique, 
analyzed with the same method, accurately and 
consistently calibrated in the same photometric 
system, and with well determined incompleteness 
factors and photometric errors. Our photometric 
sample is being built following these requirements. 

Up to now we have acquired the photometry 
for more than 30 clusters, in some cases only 
in the B and V bands, in others also in the I 
band, and in a few cases also in U and addi- 
tional bands. About half have been analyzed. All 
magnitudes have been calibrated to the Johnson- 
Cousins system, via observations of primary and 
secondary standard stars (Landolt 1992, Stetson 
2000) properly selected to cover the color range 
of the cluster stars. All the data have been re- 
duced with PSF fitting packages, particularly suit- 
able for crowded fields photometry; we have used 
DAOPHOT (Stetson 1987, Davis 1994) in most 
cases, ROMAFOT (Buonanno et al. 1983) and 
PSFex (E. Bertin, priv. communication) in a few 
others. In all cases incompleteness factors and 
photometric errors have been determined with the 
help of extensive artificial star tests. All the de- 
tails can be found in the original papers quoted in 
Table 1. 

We have already published the results for 16 
of these clusters^, which are listed in Table 1, 
together with the derived age, distance modu- 
lus, reddening, indicative metallicity and refer- 
ence. For sake of homogeneity, all the listed quan- 
tities refer to the best cases obtained with syn- 
thetic CMDs (see next Sections) based on the 
Padova 1994 stellar models (Bressan et al. 1993, 
Fagotto et al. 1994), because these are the tracks 
which more often turned out to provide the best 
fit to the data. There are (a few) cases, however, 
where other stellar models were found to better 
reproduce the data (see next Sections and the pa- 
pers quoted in Table 1 for details). 

The metallicity indicated in Table 1 is that of 

All the published data are available in electronic form from 
the BDA (Base Dcs Amas) maintained by J.C. Mermilliod 
at http;/ /obswww.unige.ch/webda/webda.html 



the stellar models in better agreement with the 
data and cannot be taken as an accurate deter- 
mination of the cluster elemental abundance, for 
which detailed spectroscopic analysis is necessary. 

The clusters position in the Milky Way is shown 
in Fig. 1, where the four main spiral arms, as delin- 
eated by Georgelin & Georgelin (1976) and Taylor 
& Cordes (1993) are also displayed. The Sun is 
assumed to be located at 8 kpc from the Galactic 
center. The; clusters already analyzed and pub- 
lished within the BOCCE photometric project are 
shown as filled dots, whilst those with ongoing 
analysis (listed in the first lines of Tabic 2) are 
represented by open circles; open triangles indi- 
cate the remaining clusters (also listed in Table 2) 
for which the photometric observations have al- 
ready been performed. 

Fig. 2 shows the empirical CMDs in the B-V 
color of the 16 clusters analyzed so far. They are 
plotted in order of decreasing age (see Table 1 and 
Section 4) rightwards from the top left panel. The 
corresponding CMDs in V-I for those clusters for 
which observations were performed also with the 
I filter are shown in Fig. 3. 

The appearance of the CMDs clearly reflects 
both the Galactic location of the cluster and the 
telescope performances. Most clusters were ob- 
served with CCDs mounted at 1 meter-class tele- 
scopes (the 1.5m Danish and the 0.9 Dutch tele- 
scopes in La Silla, Chile, the 1.5 m Bologna Obser- 
vatory one in Loiano, Italy, the Palomar 60 inches, 
USA). The exceptions can easily be recognized 
from the depth of the CMDs: NGC 6819, NGC 
2099, NGC 2168 and NGC 2323 come from the 
CFHT open cluster survey (Kalirai et al. 2001a, 
2001b, 2001c, Kalirai et al. 2003) and NGC 6939 
from the Italian National telescope (TNG) on the 
Canary Islands, both being 4 meter-class tele- 
scopes. Finally, Be 22 and Be 29, while observed 
with a similar diameter telescope (the ESO New 
Technology Telescope on La Silla), are the farthest 
ones. 

Field stars contamination appears in all CMDs. 
Its prominence depends of course on the sky 
area sampled (from the 3.8 x 3.8 arcmin^ of 
the Dutch telescope, to the 42 x 28 arcmin^ of 
the CFHT), but generally speaking, the clusters 
seen in the direction of the Galactic Anticenter 
(e.g. Be 21, Be 29 and Cr 110) are less affected 
by fore/background contamination than those lo- 
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catcd towards the center (NGC 6253 and Cr 261). 
Furthermore, the contamination is lower in the 
third quadrant (180° < 1 < 270°) than in the oth- 
ers, with nearby systems like NGC 2099, NGC 
2168, NGC 2323, NGC 6819 and NGC 6939, with 
galactic longitude between ~ 74° and ~ 221°, as 
heavily contaminated as NGC 6253 and Cr 261. 
This depends partly on the proximity of populous 
spiral arms and partly on the presence of over- 
densities such as the disk warp and/or accreted 
satellites streams (see Discussion). 

3. The application of the synthetic CMD 
method 

Age, reddening and distances of all clusters 
have been derived applying the synthetic CMD 
method (Tosi et al. 1991) to the empirical CMDs. 
The best values of the parameters are found by se- 
lecting those providing synthetic CMDs with mor- 
phology, number of stars in the various evolution- 
ary phases and luminosity functions (LFs) in bet- 
ter agreement with the empirical ones. 

The synthetic CMDs are constructed via Mon- 
teCarlo extractions of mass - age pairs, according 
to an assumed IMF, SF law, and time interval of 
the SF activity. Each extracted synthetic star is 
placed in the CMD by suitable interpolations on 
the adopted stellar evolution tracks and adopting 
the Bessel, CastelU & Pletz (1998) tables for pho- 
tometric conversion to the Johnson-Cousins pho- 
tometric system. The absolute magnitude is con- 
verted to a provisional apparent magnitude by ap- 
plying (arbitrary) reddening and distance modu- 
lus. The synthetic stars extracted for any pro- 
visional magnitude and photometric band are as- 
signed the photometric error derived for the actual 
stars of the same apparent magnitude. Then, they 
are randomly retained or rejected on the basis of 
the incompleteness factors of the actual data, de- 
rived from extensive artificial star tests. 

Once the number of objects populating the 
whole synthetic CMD (or portions of it) equals 
that of the observed one, the procedure is stopped, 
yielding the quantitative level of the SF rate con- 
sistent with the observational data for the pre- 
scribed IMF and SF law. To evaluate the good- 
ness of the model predictions, we compare them 
with: the observational LFs, the overall morphol- 
ogy of the CMD, the stellar magnitude and color 



distributions, the number of objects in particu- 
lar phases (e.g., at the main sequence turnoff, MS 
TO, on the red giant branch, RGB, in the clump, 
etc.). A model can be considered satisfactory 
only if it reproduces all the features of the em- 
pirical CMDs and LFs. Given the uncertainties 
affecting both the photometry and the theoreti- 
cal parameters (stellar evolution tracks included), 
the method cannot provide strictly unique results; 
however, it allows us to significantly reduce the 
range of acceptable parameters. 

In this way we derive the age, reddening, and 
distance of the cluster, and indicate the metallicity 
of the stellar evolution models in better agreement 
with the data. In principle, the latter should be 
indicative of the cluster metallicity, but this de- 
pends significantly on some of the stellar model 
assumptions, such as opacities. 

To estimate if, and how many, unresolved bi- 
nary systems could be present in the cluster, the 
synthetic CMDs are computed either assuming 
that all the cluster stars are single objects or that 
an (arbitrary) fraction of them are members of bi- 
nary systems with random mass ratio. The effect 
of unresolved binarism on the V magnitude and B- 

V color was first addressed by Maeder (1974) and 
more recently discussed by Hurley & Tout (1998). 
Here we have re-calculated the effect also on the 

V I color using the Padova isochrones (Bertelli et 
al. 1994) and included it in the construction of 
the synthetic CMDs. Fig. 4 shows the adopted 
magnitude and color variations of MS objects of 
solar metallicity as a function of the mass ratio be- 
tween the secondary and the primary component 
of the imresolved system. Again, the comparison 
of the resulting morphology and thickness of the 
synthetic MSs of each cluster with the observed 
ones allows us to derive information on its most 
likely fraction of binaries. 

The advantage of the synthetic CMD method 
over the classical isochrone fitting to derive the 
clusters parameters is manifest: with the synthetic 
CMDs we take into account both the number of 
stars predicted by the adopted stellar tracks in 
each evolutionary phase (directly related to the 
stellar lifetimes in that phase) and the effects of 
the photometry (errors, incompleteness and im- 
age blend) on the morphology and spread of the 
various CMD sequences. This clearly provides a 
larger number of constraints and leads to a better 
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selection of the solution, if not to its uniqueness. 
For instance, an isochrone may appear to well fit 
the data, but then the corresponding distribution 
of the number of stars on the MS or at its TO is 
completely at odds with the empirical ones and the 
case must be rejected, or be given a lower ranking 
with respect to other models (see e.g. Kalirai & 
Tosi 2004) . Moreover the blue and red edges of the 
MS distribution are important tools to quantify 
both the reddening (its mean value in the cluster 
direction and the existence of differential varia- 
tions) and the presence of binary stars. 

To test the effect of different input physics on 
the derived parameters, we have run the simula- 
tions of each cluster CMD with several different 
types of stellar evolutionary tracks. The adopted 
sets were chosen because they assume different 
prescriptions for the treatment of convection and 
range from no overshooting to rather high over- 
shooting from convective regions. They are thus 
suitable to evaluate the intrinsic uncertainties still 
related to stellar evolution models. The three 
types of stellar models applied to all clusters are: 
the tracks of Ventura et al. (1998 and private 
communication, hereafter FST) with high (rj = 
0.03), moderate (ry = 0.02) and no (77 = 0.0) over- 
shooting; the tracks of the Padova group (Bressan 
et al 1993, Fagotto et al. 1994, hereafter BBC) 
with overshooting; and the tracks of the Frascati 
group (Dominguez et al. 1999, hereafter FRA) 
with no overshooting. The parameters of some 
of the clusters were also inferred, in the origi- 
nal papers, adopting the models computed by the 
Geneva group (first Maeder & Meynet 1989 and 
then Schaller et al. 1992 and Charbonnel et al. 
1993) and by Vandenbergh (1985). 

Some of the clusters of our sample were ana- 
lyzed and published before the advent of tli(^ stellar 
models currently available and their parameters 
were then derived with the same method but with 
different stellar tracks than clusters more recently 
studied. For such systems we have rc-derivcd age, 
distance and reddening with the version of the 
BBC, FRA and FST models quoted above, so that 
all the results presented here have been obtained 
consistently. For sake of homogeneity, also the fore 
coming analyses of the remaining clusters of the 
BOCCE project will be performed with the same 
sets of evolution tracks. 

For those cases where the metallicities are not 



well constrained, we allow different values in our 
comparisons. We consider as solar metallicity 
tracks those with Z=0.02, because they are the 
ones calibrated by their authors on the Sun, inde- 
pendently of the circumstance that nowadays the 
actual solar metallicity is suggested to be lower 
(see Asplund et al. 2004). 

All models assume that the star formation 
activity has lasted 5 Myr (i.e., approximately 
an instantaneous burst, relative to the age of 
most of our studied clusters) and that the stars 
were formed following a single slope (x = 1.35) 
Salpeter's IMF over the whole mass range covered 
by the adopted tracks (0.6 - 100 Mq): simplis- 
tic assumptions which however don't affect our 
results, as demonstrated by several checks made 
assuming both different burst durations and IMFs. 

4. Ages, reddenings and distances from 
the CMDs 

Table 3 lists age (in Gyr), intrinsic distance 
modulus and reddening derived for each cluster 
with the synthetic CMD method using the three 
different sets of stellar models described above. 
Each triplet of parameters refers to the best cases 
obtained with the BBC, FRA and FST models, 
for all the tested metallicities. The differences in 
the derived quantities provide a quantitative es- 
timate of the intrinsic error attributable to the 
stellar evolution tracks. 

Age is the parameter more affected by the the- 
oretical uncertainties. Table 3 shows that there 
can be a factor up to 1.8 between ages derived 
with stellar models with or without overshooting 
from convective regions. It is well known that, 
for a given mass, stellar models taking overshoot- 
ing into account have larger cores and are there- 
fore brighter than models without it. This makes 
them appear as stars of slightly larger mass. Cor- 
respondingly, an object of given luminosity, in 
the framework of overshooting models has a lower 
mass (and therefore an older age) than with no 
overshooting models. Indeed, this is what we find 
for the majority of our studied clusters: the BBC 
and FST models provide ages older than the FRA 
models. The few exceptions to this rule are due to 
uncertainties on the best fitting model metallicity. 
If we do not consider the two oldest clusters Cr 261 
and NGC 2243, for which the BBC models appear 
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to have problems (sec next Section), the ages de- 
rived with the two types of overshooting models, 
BBC and FST, are in fairly good agreement with 
each other, since they differ at most by 20%. 

Distance modulus and reddening appear to be 
much less affected by differences in the adopted 
stellar models, specially if we take into account 
that the available BBC, FRA and FST evolution- 
ary tracks do not have the same nominal metal- 
licitios. In the range of interest for open clusters, 
the BBC models are available for metallicity Z = 
0.004, 0.008, 0.02 and 0.05, the FRA models for Z 
= 0.001, 0.006, 0.01 and 0.02, and the FST mod- 
els for Z = 0.006, 0.01, 0.02 and 0.03. Hence, 
only the CMDs of solar metallicity clusters could 
be compared with synthetic CMDs with the same 
metallicity for all types of tracks. This leads to 
an uncertainty in the derived reddening, which in- 
evitably implies an uiic;ertainty also in the distance 
modulus. On the other hand, all theoretical CMDs 
were transformed to the observational plane with 
the same photometric conversion tables (Bessel et 
al. 1998), instead of directly taking values given 
by each group, thus avoiding a further source of 
uncertainty. 

Since the BBC models have turned out to bet- 
ter reproduce the overall features of the empiri- 
cal CMDs of most of the clusters, we will adopt 
their resulting parameters as reference values, of 
age, intrinsic distance modulus (m-M)o and red- 
dening E(B— V) attributed to each cluster by the 
BBC stellar models from the application of the 
synthetic CMD method. The corresponding TO 
mass, model metallicity, distance d from the Sun 
and R from the Galactic center are also given. In 
a few cases (see next Section) the BBC evolution 
tracks were found to reproduce the clusters ob- 
served properties less satisfactorily than the FST 
or the FRA models. We have flagged with a colon 
the values in Table 1 corresponding to these cases, 
which are to be considered more uncertain. In the 
case of Cr 261, where the BBC models are likely 
to have failed. Table 1 shows the values inferred 
with the FST models. These values happen to 
be in agreement with those derived for Cr 261 by 
Carraro, Girardi & Chiosi (1999) via isochrone fit- 
ting with the new Padova stellar models (Girardi 
2000). 

Detailed comparison with other analyses of each 
cluster has been performed in the individual pa- 



pers, and we do not repeat it here; we only men- 
tion the recent work on NGC 2243 by Anthony- 
Twarog, Atwell & Twarog (2005), since our orig- 
inal paper is old, and this study constitutes an 
important improvement on this cluster. As part 
of a larger project involving a series of key open 
and globular clusters, they used Stromgren pho- 
tometry to derive the cluster parameters, argu- 
ing that intermediate-band photometry is able to 
reach more precise results than broad-band pho- 
tometry, in particular for the metallicity. NGC 
2243 was targeted because it is old, metal-poor, 
and just outside their proposed discontinuity near 
R-GC = 10 kpc. They compare their results with 
Bonifazi et al.'s (1990), finding good agreement. 
This is still true with our present analysis, since 
they find as best values E(B-V) = 0.055, [Fe/H] = 
—0.57, age = 3.8 Gyr, and apparent distance mod- 
ulus 13.15 (see our updated parameters in Tables 
1 and 3). They also rightly caution that, since the 
resulting estimates are tied to the adopted set of 
evolutionary tracks, it is safer to work on a rela- 
tive scale to derive the properties of a sample of 
clusters. They also compare Be 29 to NGC 2243, 
again finding very good agreement with our deter- 
minations. 

5. The cluster CMDs as stellcir evolution 

tests 

When the BOCCE project started many years 
ago (Romeo et al. 1989), the stellar evolution 
models available in the literature were still signif- 
icantly affected by uncertainties in the treatment 
of the convective regions, in the adopted opacity 
tables, in the normalization to observed systems. 
The photometric conversions from the theoretical 
effective temperature - luminosity plane to the 
empirical color - magnitude plane were very un- 
certain too and often created ad hoc for each grid 
of stellar tracks, to let it reproduce the colors of 
observed objects. Moreover, most stellar models 
were computed for relatively small mass ranges, 
and it was therefore quite difficult to set up a ho- 
mogeneous set of models covering the whole mass 
range of interest for open clusters (i.e. low and 
intermediate mass stars from 0.6 to 8-10 Mq). As 
a c:ons(Kiucncc, several of the adopted models were 
unable to properly reproduce all the CMD features 
of our observed clusters; first of all, the color of 
the red giant branch (RGB) and the morphology 
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of the MS TO. 

The three types of stellar models (BBC, FRA 
and FST) adopted now are those that have proven 
most successful so far in reproducing the observa- 
tional properties of the clusters. Nonetheless, none 
of them reproduces perfectly all the CMD features. 
The differences between prcidictcd and observed 
features sometimes may be attributed to uncer- 
tainties in the data^, but in most cases they are 
the systematic consequence of the stellar model 
assumptions. 

The FST models are very instructive in this 
respect because, for any given mass and metal- 
licity, they provide three sets of tracks computed 
with the same assumptions except the amount of 
overshooting. By comparing the synthetic CMDs 
based on these three different amounts of over- 
shooting (?7=0.0, 0.02 and 0.03) with the empiri- 
cal CMDs of our clusters we have found (Ventura 
et al. in preparation) that the models with over- 
shooting are in much better agreement with the 
data than the models without it. Even for Cr 261, 
where core overshooting is in practice not effective 
because at that old age (6 Gyr with the FST mod- 
els) MS stars have low mass and hence radiative 
cores, the 77=0.02 models are slightly better than 
the r]=0.0 ones. The only exception to this rule is 
NGC 2243 (which is also rather old, with an age of 
4 Gyr with the BBC models and 2.9 with the FST 
ones), where all the FST models fail to provide a 
good fit to the empirical CMD, but the ones with 
77=0.0 are less inconsistent than the others. In all 
the other cases with TO mass large enough to let 
the stellar core be convective (and overshooting 
be cfFectivo), the 77=0.0 models cither do not pop- 
ulate correctly the clump and/or the TO region, 
or do not show the correct MS curvature. 

It is interesting that, while the need of over- 
shooting in the FST models is evident, it is al- 
most impossible to systematically choose between 
7] = 0.02 or 0.03 from our data. For some clus- 

^for instance, Twcirog, Anthony- Twarog & De Lee 2003 have 
suggested the existence of a color term in the calibration of 
our photometry of NGC 6253, comparing it to other analy- 
ses of the same cluster; for this cluster there also appear to 
be quite large zero point shifts among available photome- 
tries. Anyway, the values they derive for this cluster are 
in reasonable agreement with our findings: metallicity def- 
initely higher than solar, age between 2.5 and 3.5 Gyr, red- 
dening 0.26 and apparent distance modulus between 11.6 
and 12.2. 



tors wc don't find significant differences in the 
agreement between empirical and synthetic CMDs 
switching from 0.02 to 0.03; for other clusters (Be 
29, Cr 110, NGC 6939, Pismis 2) we find that the 
77 = 0.02 models better reproduce the data; but for 
yet other clusters (NGC 6819, NGC 2660, NGC 
2168 and, possibly, NGC 2323) wc find instead a 
better fit with the tj = 0.03 models. Looking at 
the values listed in Table 1, one might be tempted 
to find hints for an anti-correlation between the 
age of the cluster and the amount of overshooting, 
but we think that the sample is still too small, the 
uncertainties too large and the relation too loose 
to support this suggestion. 

One intriguing aspect of stellar evolution mod- 
els concerns the slope of the lower MS. Sometimes 
the photometry is not accurate enough to allow 
for a safe definition of the MS shape and curvature 
at faint magnittides, but there arc cases in which 
the photometry is reliable and the synthetic CMDs 
(and isochrones) do not reproduce them perfectly. 
This problem has been recently addressed by Gro- 
cholski & Sarajedini (2003), who presented a com- 
parative study of isochrone fits to CMDs in the vi- 
sual and infrared bands. They took five sets of the- 
oretical models and compared them in the lumi- 
nosity - temperature plane. They fit the observed 
CMDs of six open clusters, of various ages and 
metallicities, with the five sets of isochrones and 
found a variety of inconsistencies on the lower MS 
slope. Their results are somewhat affected by the 
fact that they use the isochrones as transformed 
from the theoretical plane (luminosity and temper- 
ature) to the observational one (magnitudes and 
colors) by each group of stellar modelers, i.e. fol- 
lowing different prescriptions and introducing dif- 
ferential variations on the slope itself. In our case, 
this differential (and unknown) effect is absent, 
because we transform all the theoretical models 
with the same conversion tables. Still, we find the 
same kind of problems as described by Grocholski 
& Sarajedini: when the observed CMDs are deep 
enough, the fit to the lower main sequence is often 
(but intriguingly not always) bad. This happens in 
general around My — 8 for the visual colors, and 
is attributed by Grocholski & Sarajedini to some 
ingredient still missing in the model atmospheres 
for low mass stars. Indeed, this is likely to be at 
least one of the causes. For instance, the circum- 
stance that most photometric conversions are ap- 
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plied under gray atmosphere assumptions, which 
are Ukely not to be applicable to stars with My > 
8 is a probable source of inconsistency. Whatever 
the reason, there is no single set of tracks able to 
match all CMDs in all colors, at all ages, and for 
all metallicities. All of them behave quite well in 
some cases, but not in others. 

With our sample of OCs, we have found that 
the FST models in some cases predict too steep 
MS slopes at the lower masses, in spite of being 
usually able to reproduce the observed features of 
the upper MS and TO better than other models. 
In fact, we have found that, while reproducing the 
colors and magnitudes of the upper MS, the faint 
end of the FST synthetic MS is too blue (and bluer 
than that of the BBC and FRA models) for NGC 
2243, NGC 6819 and NGC 6939. In a few cases 
(e.g. NGC 2660) the FST models overpredict the 
number of faint MS objects more than other stellar 
tracks (i.e. more than attributable to evaporation 
effects). 

The FRA models are the only ones able to well 
reproduce the MS curvature of NGC 6819 and Pis- 
mis 2, and fit perfectly the older cluster Cr 261. 
On the other hand, they often overpopulate the 
post-MS phases and present a too hooked TO mor- 
phology at young and intermediate ages (see e.g. 
NGC 2099, NGC 2168, NGC 2323 and NGC 6819 
in Kalirai & Tosi 2004), which make them less vi- 
able than other tracks for clusters younger than 
~2 Gyr. 

The BBC tracks turn out to be those in over- 
all better agreement with the observed CMDs. In 
other words, they steadily provide a good fit to 
the data of (almost) any age and mctallicity, even 
if they may appear to reproduce detailed features, 
such as MS gaps and bumps, TO shape, clump 
numbers, in a less strikingly good manner than 
other models. The only case where the BBC tracks 
have not led to a convincing result is Cr 261. For 
that cluster we (Gozzoli et al. 1996) foimd an 
age of 11 Gyr with the BBC models, much older 
than with tracks without overshooting, as if over- 
shooting were strongly active in the cluster stars, 
despite their low mass (the TO mass is about 1.1 
Mq.) which implies a radiative and not a convec- 
tive core. We have reexamined its case, with the 
FRA and FST models, not available at the time 
of that study, and found a very good fit to the 
observational CMD for an age of 5 Gyr with the 



FRA models or 6 Gyr with the FST models (inde- 
pendently of the adopted amount of overshooting, 
as expected, since overshooting should not be ac- 
tive in low mass stars). The problem with the 
BBC tracks (and isochrones, Bertelli et al. 1994) 
is that at ages of 6 to 9 Gyr they show a hooked 
TO, which is usually appropriate for younger sys- 
tems and does not reproduce the smooth tighter 
morphology of the TO of Cr 261. At older ages, 
the BBC TOs recover the proper shape, but then 
the clump luminosity becomes excessively bright 
for Cr 261. It is interesting to note that the 
new Padova isochrones (Girardi 2000) present TO 
shapes much more similar to those of the FST and 
FRA models (and of Cr 261, see also Carraro et 
al. 1999 and their fig.l) at these ages. We thus 
feel that the BBC stellar models might have had 
some problems in the TO evolutionary phases of 
low mass stars, subsequently overcome in more re- 
cent computations. 

6. Field contamination, reddening and bi- 
nary systems 

A safe derivation of cluster age, reddening and 
distance from CMDs requires not only reliable 
photometry, with accurate calibration, but also a 
careful analysis of the back/foreground contami- 
nation and as much as possible information on the 
actual cluster membership of the measured stars. 
In fact, the shape and the luminosity of the corner 
stones for the parameters derivation (upper MS, 
TO and clump) may be significantly altered by the 
presence of line of sight intruders. Unfortunately, 
studies of radial velocities or proper motions of 
the individual stars are available only for a few 
open clusters and often the only means to estimate 
cluster membership is by comparison with nearby 
blank fields. De-contamination is even more cru- 
cial if one aims at deriving also the cluster initial 
mass function or the star formation rate. We thus 
emphasize the importance of acquiring the pho- 
tometry of appropriate blank fields, even if this 
requires additional observing time and a careful 
choice of the field position. The case of the beau- 
tiful CMDs from the CFHT open cluster survey 
(Kalirai et al. 2001a, 2001b, 2001c, Kalirai et al. 
2003) is illustrative. Despite the large field of view 
of the CFH12K mosaic camera (42' x 28') Kalirai 
& Tosi (2004) have found that in the case of NGC 
2168 and NGC 2323 many cluster members are 
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still present even at the edges of the observed field: 
the lack of observations of an external field pre- 
vents any reliable estimate of the back/foreground 
contamination. 

Open clusters are affected by dynamical friction 
with the Galactic disk and by stellar evaporation, 
so that stars initially members can be found at 
rather large distances from the cluster center. An 
example is given e.g. by Be 22 (Di Fabrizio et al. 
2005, see figs. 4 and 5), for which we acquired a 
comparison field at a distance of 30 arcmin from 
the cluster center that still shows a main sequence 
similar to that of the cluster. Since this MS is 
not present in the Galactic synthetic models of 
the field disk population by Robin et al. (2003), 
its most likely interpretation is that cluster mem- 
bers are present even at such large distance. This 
makes it additionally uncertain to infer member- 
ship by statistically subtracting the objects found 
in adjacent fields. For this reason, we have gen- 
erally restricted our derivation of the cluster pa- 
rameters to the CMDs of their most central re- 
gions, where field contamination is smaller. The 
synthetic CMD method, by providing the theoret- 
ical number of stars in all evolutionary phases, al- 
lows to easily individuate the clusters where evap- 
oration is more likely to have occurred. Clusters 
whose fainter main sequences are systematically 
underpopulated with respect to those of the cor- 
responding synthetic CMDs are robust candidates 
for evaporation, once the reliability of the incom- 
pleteness factors is guaranteed by extensive and 
appropriate artificial star tests. In our sample, we 
have found that NGC 6253, NGC 6819, Be 29 and 
Cr 261 have faint MSs systematically less popu- 
lated than the corresponding ones predicted by all 
the stellar models. We thus consider these clusters 
as having suffered significant evaporation of low 
mass stars during their life. These four systems are 
among the oldest ones of the provisional BOCCE 
sample and have had more time than the others 
to loose an appreciable fraction of their stars. It is 
interesting to notice, though, that Be21 and NGC 
2243, which are also older than 2 Gyr, do not show 
evidence of evaporation. NGC 2243 is by far the 
OC more distant from the Galactic plane in our 
entire sample, and this may explain it; no simple 
explanation is available for Be 21. 

We have not generally tried to determine the 
extension of our DCs, partly because it was out- 



side the main goal of our project, partly because 
our field of view is often so small as to make 
it impossible. Only for the farthest clusters (Be 
29, Be 22 - already discussed in Di Fabrizio et 
al. 2005 - Be 21) or for the few observed with a 
larger field of view (NGC 6939, Cr 261) the ex- 
ercise if worthwhile. The radial extent of the 4 
OCs observed with the CFHT has already been 
determined (Kalirai et al. 2001b,c; 2003): NGC 
6819 has a radius of 9.5 arcmin, NGC 2099 of 
13.9 arcmin, while the determinations for NGC 
2323 (about 15 arcmin) and NGC 2168 (more than 
20 arcmin) are more uncertain, since the clusters 
seem to fill the entire field of view. We computed 
the radial star density distribution for Be 29, Be 
21, NGC 6939, and Cr 261. We show in Fig. 5 the 
projected density distribution (in annuli 10 or 20 
arcsec wide) for the 4 OCs, giving distance form 
the centers both in arcsec and in parsec. The only 
one for which a flattening of the distribution is in- 
disputable, so we can determine where the cluster 
disappears, is Be 29 (with a radius less than about 
3 arcmin, in good agreement with the value in the 
Dias et al. 2002 catalogue). The distribution for 
Be 21 is still decreasing at our field's edge, and the 
same seems to be the case for NGC 6939, where 
we chose only stars brighter than V = 22, to avoid 
problems with an excessive field stars fraction and 
incompleteness of our photometry. Finally, Cr 261 
(where we considered only stars brighter than V 
= 20) seems to constitute a small fraction of the 
observed field, as already found by Gozzoli et al. 
(1996). 

As many disk objects, our clusters are gener- 
ally affected by significant reddening. However, 
only few of them clearly suffer also from differen- 
tial reddening. The synthetic CMD method, by 
distributing the stars in the CMD according not 
only to the evolutionary phase, but also to pho- 
tometric errors and fractions of presumed binary 
systems, allows one to quantify differential redden- 
ing more safely than with simple isochrone fitting. 
A varying reddening widens the whole MS by a 
given amount AE(B— V) , while unresolved bina- 
ries produce a secondary sequence redder/brighter 
than the single stars MS, intersecting it above the 
TO region. In practice, differential reddening and 
binarism are distinguishable from each other be- 
cause the former implies a smooth spread of the 
MS from the TO down to its faintest portions. 
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whilst unresolved binaries, even if distributed with 
random mass ratios, imply a split of the MS and 
the appearance of a bright appendix (up to 0.75 
mag brighter, corresponding to equal mass bina- 
ries) on top of the TO. Based on these arguments, 
we have found that Be 21 and Pismis 2 are the only 
systems of our current sample clearly affected by 
differential reddening (0.04 < AE(B-V) < 0.1). 
Work in progress (Bragaglia et al. 2005, in prepa- 
ration) on the photometric data and the spectra 
of three stars in NGC 3960 indicate also for this 
cluster some differential reddening, confirming the 
finding by Prisinzano et al. (2004). 

All the clusters, instead, seem to contain unre- 
solved binary systems, with estimated percentages 
varying between 15% and 60%. 

NGC 2099, NGC 2168, NGC 2323 and NGC 
6819, in addition to the binary sequence, show 
MS widths larger than attributable to the mea- 
sured photometric errors (Kalirai & Tosi 2004). 
This could be due to a differential reddening of 
the order of 0.01 mag and/or to internal metal- 
licity spread. However, in these cases we rather 
tend to ascribe the MS width to the difficulty in 
calibrating the different CFH12K CCD chips (Mc- 
Cracken, private communication). 

7. Discussion and preliminary results 

With the derived reddening and distance mod- 
ulus, the observational CMDs can be plotted as 
absolute magnitude vs intrinsic color (Fig. 6). 
In this reference frame, the TO brightness ob- 
viously anticorrclatcs with the cluster age, while 
the clump luminosity at first sight is roughly con- 
stant, although actually inversely proportional to 
the chemical abundance. At first order approx- 
imation, the magnitude difference, AV, between 
TO and clump is thus an age indicator. A safe 
use of this indicator requires a careful and homo- 
geneous definition of the reference TO and clump 
points, which are affected by photometric uncer- 
tainties, including blending of unresolved stars, 
small number statistics, and back/foreground con- 
tamination. In practice, the identification of the 
reference points differ from one author to the 
other. For instance, Twarog & Anthony-Twarog 
(1989) define the morphological age ratio (MAR), 
reddening independent and almost metallicity in- 
dependent, as the magnitude difference between 



the brightest TO point and the clump, divided by 
the color difference between the RGB color at the 
clump luminosity and the bluest TO point. Janes 
& Phelps (1994) adopt as morphological age in- 
dex (MAI) the magnitude distance from the red 
clump of the inflection point between the TO and 
the base of the RGB, as did more recently Salaris, 
Weiss & Percival (2004). Carraro & Chiosi (1994) 
assume that the reference TO luminosity is 0.25 
mag fainter than observed in the CMD, because of 
the possible presence of unresolved binaries which 
brighten the apparent TO. 

The oldest cluster of our provisional sample is 
Cr 261, whose age has been conservatively esti- 
mated with the FST models, due to the problems 
with the BBC tracks discussed in the previous Sec- 
tion. Its AV is 2.6. The next older systems are in 
the interval 4 - 3 Gyr (Be 29, NGC 2243 and NGC 
6253), with 2 < AV < 2.5. The clusters younger 
than 1 Gyr (NGC 7790, NGC 2323, NGC 2168, 
NGC 2099 and NGC 2660) all have AV between 
0.2 and 0.5 and both the TO and the clump are 
much more difficult to individuate, because of the 
small number of cluster members in these evolu- 
tionary phases. 

Fig. 7(a) shows that AV seems well correlated 
with the cluster age derived with our method, and 
the line represents a preliminary calibration of this 
relation. Although it appears as a very reasonable 
approximation, we will not make use of it in its 
present form: these 16 clusters are still too few to 
give a reliable calibration between AV and age be- 
cause they were not chosen for this purpose, and 
do not sample the entire relevant age range (see 
e.g. Salaris et al. 2004 for a discussion on their 
choice of calibrating clusters). In particular, we 
miss the oldest OCs, and a definitive calibration 
will have to wait until at least Be 32 (age ~ 6 
Gyr), Be 17 and NGC 6791 (both with age ~ 10 
Gyr) will be properly analyzed (for the first two 
clusters work is already in progress). Furthermore, 
precise metallicitics arc, or will soon be, available 
for (most of) the clusters in our sample, either 
from our work (e.g., Cr 261: Carretta et al. 2005; 
Be 17, Be 32, and NGC 6791 for which we have 
acquired high resolution spectra) or from the lit- 
erature (e.g.. Be 29: Carraro et al. 2004; Be 17: 
Friel, Jacobson, & Pilachowsky 2005). This is im- 
portant to test the metallicity dependence of the 
AV - age relation. As explained e.g., in Salaris 
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et al. (2004), isochroncs of different mctallicity 
have different values of AV for tlie same age, since 
the luminosities of the TO and the RC depend in 
different ways on age and metal abundance. In 
particular, lower metallicities should show a larger 
AV for the same age. In our case the situation is 
confused: we have three clusters with the same 
AV but the age difference does not seem to corre- 
late univocally with the metallicity. Larger sam- 
ples, and better defined metal abundances, have 
to be employed to decide whether and how this 
influences the calibration of this useful relation. 

Prom the derived distance moduli and redden- 
ings we calculated the Galactic position of the 
clusters, used to draw Fig. 1. NGC 6253 is the 
system closer to the Galactic center and Be 29 
the farthest one. To better visualize the rela- 
tions (or lack thereof) between age, Galactocen- 
tric distance and metallicity, we have plotted these 
quantities in Fig. 7(b,c,d). We recall, however, 
that the metallicity derived here is that of the 
BBC models in better agreement with the cluster 
CMDs and should not be taken as a precise mea- 
sure of the actual cluster metallicity. To safely de- 
termine the cluster chemical abundances accurate 
spectroscopic studies are necessary. Bearing this 
caveat in mind. Fig. 7 shows that the only relation 
clearly existing between the derived quantities is 
the mctallicity dependence on Galactocentric dis- 
tance [Fig. 7(d)], while there is no apparent re- 
lation between age and metallicity. If we don't 
take the outermost point into account, we cannot 
exclude a dependence of the cluster age on Galac- 
tocentric distance; however, the spread in the sec- 
ond panel from top is quite large and we need a 
much larger sample to better pin point this rela- 
tion (or lack thereof). Similar results were already 
found and discussed by Janes & Phelps (1994) and 
Phelps, Janes & Montgomery (1994), with a much 
larger sample of clusters, but more uncertain pa- 
rameter derivation. Carraro et al. (1998) argued 
that, after correcting the cluster metallicities for 
the radial gradient effect, their clusters do show 
an age-metallicity relation similar to that inferred 
from field stars in the solar neighborhood. How- 
ever, their corrected age-metallicity distribution 
(their fig. 8) looks as scattered as ours [Fig. 7(c)], 
with no apparent trend. Yong et al. (2005) find 
no correlation between age and metallicity in their 
compiled sample of OCs. 



Do we really see a radial abundance gradient 
in the Galactic disk? Fig. 8 shows our data^ 
compared to the homogeneous samples by Priel et 
al. (2002) and Twarog et al. (1997). As aheady 
mentioned in the Introduction, the latter, on the 
basis of their re-homogenization of OC metallici- 
ties, presented the hypothesis that the true orig- 
inal metallicity distribution in the Galactic disk 
is not described by a simple linear gradient, but 
is a step distribution, with clusters inside Rqc = 
10 kpc with solar metallicity (and a dispersion of 
±0.1 dex) and those outside with a metallicity 
0.3 dex lower, with a similar dispersion (see Fig. 
8, lower panel, which uses data from Twarog et 
al. 1997). They discussed at length the different 
metallicity distributions of clusters and field stars 
and suggested that the differences are mostly due 
to the much larger orbital diffusions of individual 
stars. Corder and Twarog (2001) modeled the evo- 
lution both of a linear and a step distribution in 
metallicity for field stars and open clusters, and 
found that a step distribution does appear as a 
linear gradient if orbital diffusion is the dominant 
factor, like for field stars. For open clusters, in- 
stead, evaporation dominates, allowing to preserve 
the original discontinuity. The large sample in 
Twarog et al. (1997) was dominated by youngish 
clusters (age < 2 Gyr) and the original step distri- 
bution (which still survives) could be recognized. 
Corder & Twarog (2001) attribute the linear gra- 
dient found e.g. by Friel (1995) to the too small 
radial baseline, with most of clusters near the dis- 
continuity (in the range 8-12 kpc). They suggest 
to concentrate efforts on objects much nearer to, 
and farther from, the Galactic center to avoid this 
effect. 

Our clusters cannot prove or disprove this de- 
scription yet: we still have too few objec;ts at any 
given Rgc, we are biased towards large ages, and 
our homogeneous metallicities come from the evo- 
lutionary tracks. Our inner clusters have predomi- 
nantly solar metallicity, and of the three OCs def- 
initely outside the 10 kpc limit (remember that 

^We have translated our photometrically determined Z val- 
ues to the corresponding [Fe/H] values, for more immedi- 
ate comparison with the two other samples. We also plot 
here the available [Fe/H]'s based on high-resolution spec- 
troscopy, that are not (yet) on a thoroughly homogeneous 
scale, since three of them come from literature (see Table 
4). 
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wc arc using Rgc.o = 8 kpc), two have [Fe/H] 
~ -0.4 (from spectroscopic analyses) and one has 
Z=Z0 (Be 22, from the synthetic CMDs method; 
sec however Villanova et al. 2005). More clus- 
ters at these Galactocentric distances should be 
studied in detail and a more precise abundance 
determination has to be performed. We list in Ta- 
ble 4 the presently available metallicities based on 
high resolution studies; 5 were determined by our 
group and 3 come from literature, so they are not 
yet on a homogeneous scale. 

In any case, a simple single-slope gradient can- 
not be the best answer: our present sample con- 
tains the very distant cluster Be 29, which has 
abundances higher than what is given by a linear 
extrapolation of the gradient found for OCs nearer 
to the Galactic center (see Fig. 8, upper and mid- 
dle panels). The same is true for Saurer A with 
Rgc ^ 19 kpc and [Fe/H] ^ -0.38 (Frinchaboy 
& Phelps 2002, Carraro et al. 2004). The fact 
that inner and outer disk clusters do not seem to 
follow the same radial abundance distribution was 
already noticed by Carraro et al. (2004), who pro- 
posed two alternatives: a different chemical evolu- 
tion history for the outer disk (e.g. Chiappini et 
al. 2001) or an extragalactic origin of Be 29 and 
Saurer A*. Yong et al. (2005), combining liter- 
ature data with high-resolution spectral analysis 
of M67 and five outer disk OCs (Be 20, Be 21, 
Be 29, Be 31, and NGC 2141), also argue that 
open clusters show a linear metallicity gradient in 
the regions within 10-12 kpc from the Galactic 
center, but have a roughly constant metal content 
from there outwards. Their favored interpretation 
is that the outer clusters are the result of star for- 
mation triggered by a major merger event. Di- 
rect capture of alien clusters from accreted dwarf 
galaxies seems unlikely, since the relatively high 
[a/Fe] ratios measured in their five clusters are not 
consistent with those measured in nearby dwarf 
galaxies (Verm et al. 2004 and references therein). 

Crane et al. (2003) and Frinchaboy et al. 
(2004) connected several open and globular clus- 
ters to the Galactic Anticenter Stellar Structure 
(CASS; e.g., Newberg et al. 2002, also called Ring, 
or Monoceros Ring). Bellazzini et al. (2004) pro- 
posed a few open clusters to be linked to the Canis 

^They associated Be 29 to the Monoceros Ring, but not 
Saurer A, at variance with Frinchaboy et al. (2004, 2005) 



Major dwarf (CMa) that may (Martin et al. 2004), 
or may not (Momany et al. 2004) be a real ac- 
creted satellite, and that may perhaps be (Martin 
et al. 2005), or not be (Pefiarrubia et al. 2005) the 
cause of the GASS. The two groups have different 
approaches. Crane et al. (2003) and Frinchaboy et 
al. (2004, 2005) use position and radial velocity of 
galactic clusters to see if they lie along the GASS, 
while Bellazzini et al. (2004) look at the CMDs of 
clusters at about the same position and distance 
of the CMa system, and compare the cluster and 
surrounding field population with the CMa one. 

A few clusters in our sample have been consid- 
ered to be associated to CMa or to the GASS, on 
the basis of their CMD and/or radial velocity. For 
To 2 (Frinchaboy et al. 2004, 2005; Bellazzini et 
al. 2004) and Be 29 (Frinchaboy et al. 2004, 2005) 
there is a strong positive indication. Be 22, on the 
contrary, has been discarded both by our work (Di 
Fabrizio et al. 2005) on the basis of a comparison 
between the CMDs of the central and control fields 
and that of the CMa overdensity, and by Frinch- 
aboy et al. (2005) because its radial velocity does 
not fit the GASS longitude- velocity trend (this re- 
mains true also considering the slightly lower ra- 
dial velocity determined by Villanova et al. 2005). 
In the case of Be 29 we have tried to apply both 
methods to our data (Tosi et al. 2004): the galac- 
tic coordinates and radial velocity (measured by 
Bragaglia, Held & Tosi 2005, and by Carraro et 
al. 2004) put it right in the GASS (see also fig. 1 
in Frinchaboy et al. 2005), but our field of view 
is too small to give a definite answer about simi- 
larities between Be 29 and CMa or dissimilarities 
with the expected galactic disk population based 
on Robin's et al. (2003) model. 

All the other clusters of our present sample are 
too close to the Galactic center and/or have ra- 
dial velocities that exclude the association with 
the Ring (see Table 4). We have anyway checked 
that no unaccounted component was present in 
our CMDs by comparison with the ones produced 
by the Robin et al. (2003) Galactic model, in each 
direction, for the appropriate field of view, taking 
into account the magnitude limit and complete- 
ness factor of our photometry. This test - negative 
in all cases - is of course truly significative only for 
the cases where the field of view was large enough 
as to sample even a marginal contaminating pop- 
ulation. 
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The importance of understanding the true ori- 
gin of outer disk open clusters is clear. For in- 
stance, Be 29, the farthest known OC, is funda- 
mental to study the mctallicity gradient, if its ori- 
gin is fully Galactic. On the other hand, if Be 29 
is demonstrated to come from another galaxy with 
a different chemical history, the interest would be 
shifted towards understanding which signatures (if 
any) characterize the accreted populations, in or- 
der to identify and separate them in the Galac- 
tic disk. The Sagittarius dwarf appear to dis- 
play chemical peculiarities (Bonifacio et al. 2004, 
Mc William & Smecker-Hane 2005), and recently 
Sbordone et al. (2005) claimed that at least one 
star, possible member of the CMa, shows chemical 
signatures (abimdance ratios of a-elements, Cu, 
and heavy neutron capture elements like La, Ce 
and Nd) more appropriate for Local Group dwarf 
galaxies than for the Milky Way. 

Both To 2 and Be 29 have been studied with 
high resolution spectroscopy; do they show simi- 
lar signatures? In To 2, Brown et al. (1996) mea- 
sured a-elements and La. For the latter they find 
[La/Fe] = 0.5; for the former [a/Fe] = +0.07 dex, 
a value in perfect agreement with many other old 
OCs (see e.g. Friel et al. 2003) but unlikely for 
dwarf galaxies, where it tends to be under solar 
(Venn et al. 2004). Carraro et al. (2004) did not 
measure heavy elements in their two Be 29 targets, 
but have a-elements, that are said to be enhanced 
with respect to the solar values, at variance e.g. 
with NGC 2243, which has similar overall metallic- 
ity ([Fe/H] around —0.4 dex). However, from the 
values published in their tab. 5, [a/Fe] is about 
0.08 dex, exactly the average found for NGC 2243 
by Gratton & Contarini (1994), and again a com- 
mon value for old open clusters. 

8. Summary 

To summarize the mid-term results of the pho- 
tometric BOCCE project: 

• We confirm the tight relation between clus- 
ter age and magnitude difference between 
MS TO and clump luminosity levels, but 
caution against its application to young ob- 
jects, where small number statistics and field 
contamination may make the determination 
of the two luminosity levels too uncertain. 



• We confirm the anticorrclation between 
metallicity and distance from the Galactic 
center. We don't find evidence of discontinu- 
ities around 10-12 kpc from the center, but 
our sample is still too small to distinguish 
between a smooth gradient and alternative 
distributions. 

• We find definitely no correlation between 
cluster age and metallicity. 

• We cannot completely exclude an age- 
Galactocentric distance relation in the inner 
10-15 kpc. 

Most of the oldest clusters of the sample (NGC 
6253, NGC 6819, Be 29 and Cr 261) show clear 
evidence of evaporation of their lowest mass stars. 
Almost all our clusters appear to have stars (still 
bound members ?) with the same CMD morphol- 
ogy at quite large cluster radii. This circumstance 
makes it crucial to always observe control fields at 
appropriate distances from the cluster center. 

We emphasize the importance of deriving the 
cluster parameters with more than one set of stel- 
lar evolution models and with the same photomet- 
ric conversion tables to allow for a direct compar- 
ison of the different model predictions and for a 
correct evaluation of the theoretical uncertainties 
on the results. We have not found a set of evo- 
lutionary tracks able to reproduce in detail all 
the observed features of old and young, metal 
rich and metal poor clusters. However, all the 
three evolution sets adopted in this project per- 
form quite well in the comparison with the obser- 
vational CMDs and LFs of the examinc^d clusters, 
with the BBC models being often found in better 
agreement with the data. Overshooting from con- 
vectivc regions appears to provide a better fit to 
the observed properties of stars in the various evo- 
lutionary phases, with some hints that it might be 
more conspicuous in younger systems. As found 
also by Grocholski & Sarajedini (2003), the CMD 
region more difficult to reproduce is the lower MS, 
where stellar models often fail to provide the ap- 
propriate shape or slope. We agree with them 
that the most likely explanation for this problem 
is the inadequacy of the adopted photometric con- 
versions, based on gray atmospheric models, for 
the lowest mass stars. 
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These results are by no means conclusive, since 
we need to at least double our clusters sample 
to get enough statistics in each age and distance 
bin. Moreover, wc need to derive accurate chemi- 
cal abundances from high resolution spectra for a 
through interpretation of the apparent metallicity 
distributions. Nonetheless, we consider it useful 
to illustrate the importance of a homogeneous ap- 
proach in the derivation of the global cluster prop- 
erties and emphasize the need of further accurate 
studies of these Galactic objects which can tell us 
so much more than others on the Milky Way evo- 
lution. 
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Table 1 

Cluster parameters derived from the BBC models. We assumed Rgc,© = 8 kpc. 



Cluster 


1° 


b° 


age 


(m-M)o 


E(B-V) 


rtiTo 


Z 


d 


R 


AV 


ref 








(Gyr) 






(Mt?) 




(kpc) 


(kpc) 






Cr 261^''' 


301.68 


-5.63 


6.00: 


12.20: 


0.30: 


1.10: 


0.020 


2.74 


6.96 


2.60 


3 


NGC 2243'' 


239.50 


-17.97 


4.00: 


12.80 


0.06 


1.20: 


0.006: 


3.45 


10.20 


2.20 


2 


Be 29 


197.98 


8.03 


3.70 


15.60 


0.12 


1.20 


0.004 


13.05 


20.81 


2.00 


12 


NGC 6253 


335.46 


-6.25 


3.00 


11.00 


0.23 


1.40 


0.050 


1.58 


6.60 


2.00 


4 


Be 22 


199.88 


-8.08 


2.40 


13.80 


0.64 


1.50 


0.020 


5.75 


14.02 


1.50 


13 


Be 21 


186.84 


-2.51 


2.20 


13.50 


0.78 


1.40 


0.004 


5.01 


12.99 


1.80 


6 


NGC 6819 


73.98 


+8.48 


2.00 


12.20 


0.12 


1.60 


0.020 


2.72 


7.71 


2.00 


11 


NGC 2506 


230.56 


+9.94 


1.70 


12.60 


0.05 


1.70 


0.020 


3.26 


10.38 


1.60 


5 


Cr 110 


209.65 


-1.98 


1.70 


11.45 


0.57 


1.60 


0.004 


1.95 


9.74 


1.65 


9 


NGC 6939 


95.90 


+12.30 


1.30 


11.30 


0.34 


1.85 


0.020 


1.78 


8.37 


1.55 


10 


Pismis 2 


258.85 


-3.34 


1.10 


12.70 


1.29 


1.95 


0.020 


3.46 


9.31 


0.85 


8 


NGC 2660 


265.93 


-3.01 


0.95 


12.30 


0.40 


2.10 


0.020 


2.88 


8.69 


0.50 


7 


NGC 2099 


177.64 


3.09 


0.43 


10.50 


0.36 


2.70 


0.008 


1.26 


9.26 


0.40 


11 


NGC 2168 


186.59 


2.22 


0.18 


9.80 


0.20 


4.00 


0.008 


0.91 


8.91 


0.30 


11 


NGC 2323 


221.67 


-1.33 


0.12 


10.20 


0.22 


4.50 


0.020 


1.10 


8.85 


0.50 


11 


NGC 7790 


116.59 


-1.01 


0.10 


12.65 


0.54 


5.00 


0.020 


3.39 


9.99 


0.30 


1 



''Cr 261 is the only cluster of the Table for which the listed parameters are derived with the FST models rather 
than with the BBC ones. See text for details. 

''Colons indicate more uncertain values, resulting from non satisfactory synthetic CMDs 

References. — (1) Romeo et al. 1989; (2) Bonifazi et al. 1990; (3) Gozzoli et al. 1996; (4) Bragaglia et al. 1997; 
(5) Marconi et al. 1997: (6) Tosi et al. 1998; (7) SandrcUi ct al. 1999; (8) Di Fabrizio et al. 2001; (9) Bragaglia & 
Tosi 2003; (10) Andreuzzi et al. 2004; (11) Kalirai & Tosi 2004; (12) Tosi et al. 2004; (13) Di Fabrizio et al. 2005. 
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Table 2 

FORECOMING BOCCE CLUSTERS; DATA TAKEN FROM DiAS ET AL. (2002). On THE FIRST SIX WORK IS 
IN PROGRESS, WHILE FOR THE OTHERS WE HAVE ONLY ACQUIRED THE DATA. 



Cluster 


r 




b° 




age (Gyr) 


NGC 3960 


294, 


.37 


+6, 


.18 


1.3 


Be 17 


175, 


.65 


-3, 


,65 


12.0 


Be 32 


207, 


.95 


4. 


,40 


3.4 


To 2 


232, 


.83 


-6, 


,88 


1.0 


Mel 71 


228, 


.95 


+4, 


,50 


0.2 


NGC 4815 


303, 


.62 


-2, 


,09 


0.2 


King 11 


117, 


.16 


+06, 


.48 


1.1 


NGC 1193 


146, 


.81 


-12, 


,16 


7.9 


NGC 2204 


226, 


.01 


-16, 


.11 


0.8 


NGC 2477 


253, 


.56 


-05, 


,84 


0.7 


NGC 2849 


265, 


.27 


+06, 


,36 


0.6 


NGC 6134 


334, 


.92 


-00, 


,20 


0.9 


NGC 6603 


12, 


.86 


-01, 


,31 


0.2 


Be 19 


176, 


.90 


-03, 


,59 


3.1 


Be 20 


203, 


.48 


-17, 


.37 


6.0 


Be 66 


139, 


.43 


+00, 


,22 


5.0 


IC 361 


147, 


.48 


+05, 


,70 


0.1 


IC 1311 


77. 


.69 


+04, 


.28 


0.4 


IC 4651 


340, 


.09 


-07, 


,91 


1.1 


IC 4756 


36, 


,38 


+05, 


,24 


0.5 


Trumpler 5 


202, 


.86 


+01, 


.05 


5.0 
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Table 3 

Age, distance modulus and reddening from different stellar models 



Cluster 




age (Gyr) 






(m-M)o 






E(B-V) 




BBC 


FRA 


FST 


BBC 


FRA 


FST 


BBC 


FRA 


FST 


Cr 261 


9.00 


5.00 


6.00 


12.0 


12.2 


12.2 


0.25 


0.29 


0.30 


NGC 2243 


4.00 


3.20 


2.90 


12.8 


12.8 


13.1 


0.06 


0.08 


0.04 


Be 29 


3.70 


2.30 


3.50 


15.6 


15.8 


15.8 


0.13 


0.15 


0.10 


NGC 6253 


3.00 


3.00 


3.20 


11.0 


10.8 


10.9 


0.23 


0.33 


0.32 


Be 22 


2.40 


2.10 


2.50 


13.8 


14.0 


13.9 


0.64 


0.65 


0.65 


Be 21 


2.20 


2.50 


2.50 


13.5 


13.6 


13.5 


0.78 


0.74 


0.74 


NGC 6819 


2.00 


1.80 


2.30 


12.2 


12.0 


12.0 


0.12 


0.15 


0.10 


NGC 2506 


1.70 


1.60 


1.80 


12.6 


12.6 


12.7 


0.00 


0.07 


0.04 


Cr 110 


1.70 


1.20 


1.50 


11.5 


11.7 


11.5 


0.57 


0.52 


0.56 


NGC 6939 


1.30 


1.00 


1.30 


11.3 


11.3 


11.4 


0.34 


0.38 


0.36 


Pismis 2 


1.10 


1.10 


1.20 


12.7 


12.5 


12.7 


1.29 


1.29 


1.29 


NGC 2660 


0.95 


0.75 


1.10 


12.3 


12.3 


12.1 


0.40 


0.42 


0.37 


NGC 2099 


0.43 


0.59 


0.52 


10.5 


10.5 


10.4 


0.34 


0.16 


0.36 


NGC 2168 


0.18 


0.11 


0.20 


9.80 


9.80 


9.60 


0.20 


0.22 


0.20 


NGC 2323 


0.12 


0.11 


0.13 


10.2 


10.0 


10.0 


0.22 


0.22 


0.22 


NGC 7790 


0.10 


0.05 


0.10 


12.7 


12.7 


12.7 


0.54 


0.54 


0.54 
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Table 4 

Heliocentric radial velocities and metallicities from high resolution studies for the 

PRESENT SAMPLE. REFERENCES ARE FOR BOTH RV AND [Fe/H] WHEN THE LATTER IS AVAILABLE. WE 
also give for comparison Z, the PHOTOMETRIC METALLICITY, TAKEN FROM TABLE 1; Z = 0.004, 

0.006, 0.02, 0.05 MEAN [Fe/H] 0.7, -0.5, 0.0, +0.4, respectively 



Cluster 


RV 


[Fe/H] 


Z 


Ref 


Cr 261^ 


-25.9 


-0.03 


0.020 


1 


NGC 2243* 


+61.0 


-0.48 


0.006 


2 


Be 29 


+24.6 


-0.44 


0.004 


3 


NGC 6253'^ 


-26.2 


+0.36 


0.050 


4 


Be 22 


+95.3 


-0.32 


0.020 


5 


Be 21 


+12.4 


-0.54 


0.004 


6 


NGC 6819^ 


+4.2 


+0.09 


0.020 


7 


NGC 2506'* 


+83.7 


-0.20 


0.020 


8 


Cr 110 


+40.0 






9 


NGC 6939 


-19.0 






10 


Pismis 2 


+49.8 






11 


NGC 2660 


+21.2 






12 


NGC 2099 


+7.7 






13 


NGC 2168 


-8.0 






14 


NGC 2323 


+9.0 






15 


NGC 7790 


-78.0 
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'*[Fe/H] values measured by our group 



References. 1: Carretta et al. (2005); 2: 
Gratton & Contarini (1994); 3: Carraro et al. 
(2004); 4: Carretta et al. (2000); 5: Villanova 
ct al. (2005); 6: Hill & Pasquini (1999); 7: Bra- 
gaglia et al. (2001); 8: Carretta et al. (2004); 9: 
E. Pancino, priv.comm., preliminary; 10: Milone 
(1994); 11: Fricl ct al. (2002); 12: propr.data, 
preliminary; 13: Mermilliod ct al. (1996); 14: 
Barrado y Navascues et al. (2001); 15: Harris 
(1976); 16: Gorynya et al. (1992) 
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Fig. 1. — Positions on the Galactic plane of the 
OCs in our sample; filled circles represent pub- 
lished OCs, open circles work in progress, and 
open triangles the remaining clusters. We show 
the positions of the Sun and the Galactic center, 
and a sketch of the spiral arms (respectively: 1 = 
Norma, 2 = Scutum-Crux, 3 = Sagittarius-Carina, 
4 = Perseus) 



23 



I I 1 1 I I I I I I I I I I 1 1 I I . 
Cr ESI {<e'l_ 



I I I I I I I I I I I I I I I I I I 
NGC 3243_ 



1 1 I I I I I I I I I I 1 1 I I 1 1 
Be S9. 



I I 1 1 I 1 1 I 1 1 I I 1 1 I 1 1 I J 
NGC 6353 J 



10 



V 



15 

20 h 
10 




I II I I I II I I I II Mil l 1 1 II II I I I II I I LI I I II I 1 1 II II I I l li l II 1 1 I II I II 1 1 I 



1^- 




Be 23. 



Be 3L 



NGC SBIS 



NGC 2506. 



V 



15 h 

20 

25 
10 

15 

20 

35 
10 

15 

30 1^ 
35 



: i- . >■ 



I II I I I II I I I II ni l I I II II I I I II I I 




II 1 1 l li l II 1 1 I II I II 1 1 I J 



Cr 110. 



NGC 6939. 



PiBmiB 3. 



NGC 2660. 




Fig. 2.— V, B-V CMDs of the 16 OCs, shown in 
decreasing age order. 
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Fig. 3.— V, V-I CMDs of the 10 OCs with I data, 
shown in decreasing age order. 
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Fig. 4. — Effect of unresolved binaries on the ap- 
parent magnitude and color as a function of the 
mass ratio between the secondary and primary 
component of the system and for various masses of 
the primary star. Line symbols refer to primaries 
of the following initial mass: 26 (solid), 14 (dot- 
ted), 4.5 (short-dashed), 2.55 (long-dashed), 1.76 26 
(short-dash-dotted), 1.14 (long-dash-dotted), 0.95 
Mq (short-dash- long-dashed). The shown varia- 
tions refer to solar metallicity MS stars from the 
Padova. isochrones (^Rertelli et al. 1994V 
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Fig. 5. — Projected density distribution for Be 
29, Be 21, (both computed with a bin 10 arcsec 
wide), NGC 6939, and Cr 261 (both with bin = 
20 arcsec). 
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Fig. 6. — CMDs for the 16 OCs in absolute mag- 
nitude and intrinsic color. 
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Fig. 7. — For the clusters listed in Table 1 we 

show the distributions of: (a) age vs magnitude 
difference between TO and clump (filled circles, 
open circles, filled triangle indicate Z solar, sub- 
solar or supcrsolar, respectively), with a line rep- 
resenting our preliminary relation (see text); (b) 
age vs Galactocentric distance (same symbols); (c) 
age vs model mctallicity; (d) model metallicity vs 
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Fig. 8. — Galactocentric metallicity distribution. 
Here we show [Fe/H] instead of Z, but on the 
right side of each figure we indicate the level of a 
few representative Z values for immediate compar- 
ison. Top panel: our sample, where filled squares 
are photometric metallicities, and filled dots are 
[Fe/H] based on high resolution spectroscopy while 
open squares are the corresponding photometric 
values. Middle panel: Friel et al.'s (2002) data, 
with the gradient based on their data, which has a 
slope of -0.059 dex kpc~^. Bottom panel: Twarog 
et al.'s (1997) data, that they interpret as a step 
distribution, with the two metallicity levels indi- 
cated by the dashed fines: [Fe/H] = for R < 10 
kpc, and -0.30 for R > 10 kpc. 
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